Hopping between Random Locations: Spectrum and Instanton 
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Euclidean random matrices appear in a broad class of physical problems involving disorder. The 
£h ■ problem of determining their spectra can be mapped, using the replica method, into the study of 

a scalar field theory with an interaction of the type . We apply the instanton method to study 
their spectral tails. 
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I. INTRODUCTION 

c : 

Place N points randomly in a d-dimensional Euclidean space of volume V. Denote the locations of the points by 
Xi (i = 1, N) . Choose a suitably well-behaved function / (x), vanishing as its argument tends to infinity. Consider 
^ ! the N x N matrix 

Ht^ffr-Hj) (1.1) 

and solve the eigenvalue equation £\ Hijipj — Eipi. We are interested in the density of eigenvalues p{E) and the 
localization properties of the eigenvalues ipi as we average over the ensemble of matrices [called Euclidean random 
matrices in Ref. ( w)} generated by placing Xi randomly. The limit N — > oo, V — > oo, with the density p = y (not to 
be confused with p{E) of course) held fixed, is understood. In Ref. ([!]) p{E) was calculated in various approximations. 
(A more involved version of ( |1 . 1| ) was also studied.) We will refer to this henceforth as model J-, p, 
This type of random matrix problem may be relevant to a broad class of physical situationsBB, structural glassesta 
y—( \ and amorphous semi-conductors for example.Q Matrices of this type have also, appeared in the instantaneous normal 
^ ■ mode analysis in the theory of liquidO and in the bipartite matching problemEm in combinatorial optimization. 

A venerable problem in the study of disordered systems is that of an electron in a metal moving in a disordered 
array of impurities. The simplest version of this problem ignores the periodic potential of the metal and treats the 
' impurities as 5-function scatterers with Hamiltonian: 

\Q . 

O N 

O! H = -V 2 + 2naJ2 S (x ~ Xi) 



where Xi are randomly [located. The density of states in the tail of the distribution, where it is cxponqatially small, 
was studied by Lifshitz.EI Later this calculation was reproduced and extended using an instanton methodtJ in the case 
of a repulsive potential—a > 0. A related problem of an electron moving in a white noise Gaussian random potential 
was studied by Cardy.liil One of ustHl extended this model to include a magnetic field B in which case V 2 is replaced 
by (V — ieA) 2 with A the vector potential. The density of eigenvalues in the tail of the distribution was calculated. 
We will refer to this henceforth as model II. 

In this note, we will elucidate the relation between model I and model II. It turns out that the relationship only 
holds for the case of an attractive potential, a < 0. We will apply the instanton method to study the density pf 
eigenvalues p(E) for large negative E. This calculation is similar to that for a > already carried out earlier.EHl 
Finally, we extend model I to include a magnetic field B for d = 2. Model I with a magnetic field may be relevant for 
studying the quantum Hall transitions. 

II. REVIEW OF FIELD THEORETIC FORMULATION 

In Ref. ( the problem was mapped into a quantum field theory. We will review the procedure here. 
We start with the replica identity 
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i=l a=l / 

with = f(xi,Xj). Here a = l,2,...,n is the replica index. We have used complex integration variables <f>j here 
because eventually we wa nt to study the problem with a magnetic field in which case H will be Hermitian rather than 
real symmetric as in (1.1 ), for which real integration variables would have sufficed. Once we have calculated S_/v we 
can obtain the desired Green's function by differentiation: 

<?(*) - = **~o 5*) (2.3) 

Let us insert the representation of the identity (here ^ and f/ 1 denote two complex scalar fields) 



(2.4) 



into the functional integral refining S, thus obtaining 
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^)-(nf-) e-i^^A (2.7) 



where 

J = jd(f> e r z T^J^- l llJ^'^+ h ^ (2.6) 

is a functional of ijj. 

Integrating out 4>, we find 

2tt 

Recall that the average (J) means / ^^^-...^r-J we see that the multiple integral factorizes and we obtain 

(J) = A N (2.8) 

with 

A=i Jdxe-mW)\ a (2.9) 

We have taken the n — ► limit wherever we are allowed to do so. The integral over t/j can be done immediately giving 

E N = J D4>e$ $ dxdvr{x)rl{x - v)lp[v) A N (2.10) 

To put the factor A N into the exponential we follow Gibbs and introduce the grand canonical ensemble 

oo N . -iY, i^wi 2 

E^= / iW/* ,/_Vf (2.H) 

AT=0 ' ^ 

In other words, instead of focusing on the original problem of studying a random N by N matrix we now c onsider 
an ensemble of such matrices with N varying over the non-negative integers. We expect the sum in (2.11 ) to be 
dominated by some values of N: 

(N) = ^-logZ(a) = — I [ dxe-i^J^ x)l2 \ ► a (2.12) 

da V \.l / n^Q 
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Defining the density of points as p = ^p-, we obtain Z(a) = J Dtpe s ^ with the action 



(2.13) 



The action S(ip) defines a non-local field theory. Up to this point, any f(x, y) could have been used. A particularly 
convenient choice is the Yukawa function 



f(x,y) = f(x-y) = (-) 



d d k 



Jk{x-y) 



{2ir) d k 2 + m 2 

The overall minus sign is included so that the resulting field theory would have the standard kinetic energy term: 



(2.14) 



c r n 

Sty) = J d d x E(|V^„| 2 + m 2 |V | 2 )-pe-iE: =1 l^l 2 



One remark about the spectrum: The eigenvalue equation implies that for a general f(x) — J A^h f(k) 
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(2.15) 



(2.16) 



Thus, for f(k) < 0, as is the case with the choice in ( 2.14 ), the eigenvalues are all negative. 

A few remarks about renormalizability and dimensionless parameters are in order. This field theory is ultraviolet 
finite in d = 1, renormalizable in d = 2 (although requiring an infinite number of counter-terms in contrast to the sine 
Gordon theory with its special symmetry) and non-renormalizable in d > 2. The original problem, in model I, is well 
defined for any well defined f(x). However, for the particular choice of f(x) — (— ) j k^+m' 2 we usec ^ m OT ^ eT 

for the field theory to have a normal kinetic energy term, /(0) does not exist for d > 1. One way of regularizing is to 
write 



m = {-) 



d d k 
(2^p' 



r) = H 



d d k 



(2Tr) d (k 2 + m 2 ){k 2 + M 2 ) 



(M 2 -m 2 ). 



(2.17) 



The field theory would then be cut off correspondingly by the Pauli-Villars mass M. The ultraviolet difficulties in 
model II can be understood as arising from the fact that a ^-function potential is too singular to have a well-defined 
spectrum for d > 1. Physical ultra-violet regularizations would replace the S- function potential by a smooth function. 

Thus, in d = 1 the density of states should be a well-defined function of E , p and m. In d — 2 it should be 
expressible in terms of renormalized parameters but in d > 2 it will be strongly dependent on the details of the 
cut-off. 

Ignoring cut-offs, we see that the field theory (2.15) contains two dimensionless quantities: 

v = p/m d 

and 

tt = \E\m 2 /p, 

where we have set z = —\E\ in light of an earlier remark. The parameter v has the physical interpretation of the 
number of points in the correlation volume of the function /. Consider expanding the action: 

/n n 
d^XdVVal 2 + (m 2 - ^-Ma\ 2 ) ~ -^£>a| 2 ) 2 +■••]■ (2-18) 
a=l II II o=l 

The coupling constant of the ip 2n term is of order p/E n . Since this has dimensions of (mass) d_ ™^~ 2 \ we see that the 
condition for the dimensionless coupling constant to be small is (for n > 2): 

n n v n-i >:> L ( 2 1Q ) 

Here we have assumed O > 1, a condition neccessary for the perturbative stability of the theory, as discussed in Sec. 
IV. If is only slightly greater than 1 then we get more complicated conditions: 

d/2+n(l-d/2) 

>> 1. (2.20) 



1 - 



We see that the conditions of Eq. (2.19) or ( 2.20D require at least one of the two parameters fl and v to be large. 
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III. RELATIONSHIP BETWEEN MODELS I AND II 



Consider model II in the case of an attractive potential, a < and assume that all points are far apart compared 
to the range of the ground state wave-function, ipo(r), for a single 8- function potential. Then the lowest energy 
states will be formed by tunnelling processes between these lowest bound states. (As mentioned above, for d > 1, 
the 5-function potential must actually be replaced by some smooth function.) To formalize this, we go to the tight 
binding approximation. The hopping amplitude for an electron to go from a potential well at the origin to a potential 
well a distance R away is given by the overlap integral 



where if>o(r) is the solution of 



t(R) = 2TraM0)MR) 



(3.1) 



(3.2) 



(for r > 0) with E the (negative) binding energy in the single site potential. The hopping amplitude t{R) is thus 
just a constant times a wave function ^po(R) and so essentially equal to the Yukawa function in Eq. ( [2.14 ) with \E\ 
playing the role of to 2 . 

In ref . ( |lCl(L2|) , it was shown that model II can be represented in terms of a field theory with the action 



S(sl>) = 



/ d d x ^(|W>a| 2 - E\^ a \ 2 ) + p(l - e~ 2 ™TZ=^\ 2 

\_a=l 



In light of the preceding discussion, this is precisely what we would expect. Thus, we see that model II corresponds 
to model I with a specific choice of the function /(x), with the map of the complex variable z to the inverse of the 
strength of the potential 2-kq, and m 2 to — E. 



IV. INSTANTON ANALYSIS 



In Ref. ( [l]) the densitv-pf eigenvalues was studied starting with the field theory( 2.15 ). Here as promised we will 
use the instanton methodic to study the tail of the spectrum. 

In order for G(z) to have an imaginary part the functional integral defining the field theory has to be "sick" . 
Otherwise, the functional integral, if well defined, is manifestly real. Let us check this statement. Consider the regime 
in which z is real positive so that z = \z\. Then the potential 

V(ip) = m 2 \^\ 2 - p (e^ IV ' 12 - l) ~ (to 2 + j^J iVf + . . . (4.1) 

is well-behaved at large |?/>|and so perturbation theory should be fine. Indeed, by the argument given earlier we expect 
all the eigenvalues to be negative, and so G(z) should not have an imaginary part. 
In contrast, for z negative, we have 

V(ip) = m 2 \^j\ 2 ~ p (e^ |,/>|2 - l) ~ (to 2 - -J-J |</>| 2 + . . . (4.2) 

The potential is unbounded below for large \ip\. Perturbation theory fails and G{z) could well have an imaginary 



part. From (4.2) we see that we should distinguish two regimes: small eigenvalue \z\ < p/m 2 and large eigenvalue 
\z\ > p/m 2 . In the large eigenvalue regime the potential V{ip) starts out concave upward before becoming unbounded 
and so there can be an instanton configuration connecting tp = to tp = where V(tpo) = 0. In the small eigenvalue 
\z\ < p/m 2 regime the potential V(ip) starts out downward and the instanton approach does not apply. Some other 
non-perturbative method is needed to study the spectrum. 

To summarize, we have three regimes: (I) z real positive, no eigenvalue and G(z) does not have an imaginary part, 
(Ha) z real negative and \z\ < p/m 2 (that is, < 1) non-perturbative regime, and finally (lib) z real negative and 
\z\ > p/m 2 (that is, f l > 1 ), and the instanton approach applies provided that the coupling constants are small, 
the conditions of Eq. ( [2.19 ). In this paper we focus on the regime (lib) and hence study the tail of the eigenvalue 
distribution. 
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We will now give a heuristic argument on how large negative eigenvalues can occur. Suppose that in the distribution 
of the N random points we have an isolated cluster of k points within a length scale comparable to the length scale 
I characteristic of f(x) (I = 1/rn in our specific example). The Hamiltonian H then contains a k by k block whose 
entries are of order /(0), thus giving us one eigenvalue of order fc/(0) (recall that /(0) is negative in our example.) The 

probability Pk of obtaining such a cluster is given by the Poisson distribution Pk oc e~~( pl ) ^ p \} and so heuristically 
we obtain the estimate for the probability of obtaining a large negative eigenvalue E to be 



P{E) 



s -ItwI 1o sIjwI 



(4.3) 



We expect that this heuristic argument will work best for d = 1 . 

For z = E large and negative, let us scale tp — y/\E\(p and x — y/m to rewrite the action as 



S(<p) = 



\E\ 



-.d-2 



d d y(\VyV\ 2 + (\v\ 2 - 



Em 2 



M 2 



In terms of the dimensionless quantities ft = and v = p/m d ,we see that 



e -s _ e -unh(n) 



for some function h. 

When the conditions of Eq. ( [2.19 ) are satisfied, the functional integral giving Z(a) is dominated by the extremum 
of the action S, namely the instanton. We assume that the instanton is spherically symmetric, so that q = \ip\ is a 
function of \y\ only. Following standard practice in instanton analysis we identify q as the position of a particle along 
a line and t = -^\y\ as time, we have the equation of motion 



cPq 
dt 2 



(d-1) 



dq 
It 



dU 
dq 



(4.4) 



with the potential 



U{q) 



n 



(4.5) 



as shown in figure (1). 




FIG. 1. Potential for the single particle problem for Q — 10. 

For d > 1, there is a time-dependent friction term. For d = 1 there is no friction. 
Let us define qo and qWnby U (go) = and U' (q m in) = respectively: 



q 2 ■ 

Until 



log ft 



(4.6) 



and 



5 



~ log + log log ft - 
and ending at q 



We have q 2 
q = qo at t 
trajectory. 

Clearly, we need to treat the cases d = 1 and d > 1 separately. 

For d = 1, we have explicitly 5 = 2V2|£|m J °° <ft (±<j 2 - U(q)) 
&r — / °° dt Qg 2 — U(q)). With the initial conditions specified, the equation of motion integrates to q 2 = —2U(q) 
and so the trajectory is given by 



• • for fi >> 1. We wish to determine the trajectory of the particle, starting at 
at t = oo. In the end, we would like to calculate the action associated with the 



It is convenient to work with the reduced action 



t = 



dq' 



y/-2U{q') 



Evaluated for this trajectory ( |4.7| ) the action is equal to 

dtq 2 = 



dqy/-2U(q) 



which can be easily integrated numerically as a function of f2. 

Our instanton analysis thus predicts that for d =1, the tail of the eigenvalue distribution should go like 



P(E)~LF(\E\,p)e 



-2 v / 2^ns r .(n) 



(4.7) 



(4.8) 



(4.10) 



where L denotes the size of the system (which we need to compare with direct numerical diagonalizations of H) 
and F(\E\,p) the infamous determinant factor in instanton calculations (which we have not computed.) We expect 
.F(|-E|,p) to be slowly varying compared to the exponential factor e -2v2i/ns r (n) _ 

We have done some numerical work in which we diagonalized the matrix H directly for N ranging up to 1000. In 
figure (2) we show the numerical data for N = 1000, p = 1.2, and m — (1.3) -1 . Only the tail of the di strib ution, 
which wc take to be comprised of the 250 eigenvalues with the most negative E, is displayed. We take ( 4.10 ) and 
compute the integrated number of eigenvalues N(E) = N J_ oo dE'P(E'), treating F(\E\,p) as a constant C in the 
range of E of interest. (N(E) is not to be confused with N of course; -ZV(oo) = N.) We then do a one parameter fit 
in C to the numerical data. As expected, the theoretical curve appears to fit the numerical data in the applicable 
regime, Eq. (2.19). It should be kept in mind that the theoretical curve fails to make sense as ft approaches 1 which, 
for our particular parameter choice, corresponds to E = —2.03. Also, the numerical data is clearly dominated by 
finite size effects in the extreme tail where N(E) is 0(1). 




-3 -2 -1 

FIG. 2. The number of eigenvalues with energy < E, plotted versus E. Only the low energy tail of the distribution is plotted. 
Analytically, we are able to study the problem only for Q >> 1 in which case 

S r — —J= log fi 

v 2 



(4.11) 



We can check this result by studying the trajectory of course. Think of the particle being released at go with zero 
velocity. We divide the trajectory into three regimes: (1) short time, during which the particle moves rapidly from 



G 



So to q m im (2) intermediate time, during which the particle passes through the minimum of the potential at q m im 
and (3) long time, during which the particle slowly climbs the hill to q — 0. A priori, it is not immediately obvious 
whether the short or long time regime gives the dominant contribution to 5 r . A detailed calculation, which we now 
outline, shows that the long time regime dominates for CI ^> 1. 

In the short time regime, we approximate U(q) ~ 2qg(q — qo)) thus obtaining q{t) = qo — \q^ 2 + ' ' ' ■ Similarly, 
we approximate U (q) in the other two regimes and match the trajectory, thus obtaining in the intermediate time 
regime q(t) = q m i n — A= sin(^/2 log Cl(t — t m i n )) + • • ■ and in the long time regime q(t) = qminC^ tmin ^ + ■ ■ ■ where 

i 

t m in = (^(so-gWra))^ -> - lo f g 2 ■ Substituting into S r , we find that S r goes like 0((loglogSl)i), 0((loglogfi)^), 
and -^logCl in the short, intermediate, and long time regime respectively. Thus, the long time regime dominates 



S ~ 2\E\m\og^^- (4.12) 

P 



the heuristic argument of Eq. ( O ) 



Using /(0) = — l/2m, in d = 1, fro m Eq. (2.14), we see that the instanton result, Eq. (4.12), is in agreement with 



The limit fi = 1 + S, with 5 <C 1 is also interesting. We have qo ~ \/26 and q m in — VS, with U(q) ~ — Sq 2 + \q A 



The trajectory is explicitly found to be q(t) — y25(l — tanh y25t). We obtain 

We now turn to higher dimensions d = 2 or 3. The equation of motion now includes a friction term 

£q (d-l) dq _ 1 , 

dt 2 + t dt q{ n ' 

At small time, near qo we have 4jf + ^ ~ —7 and the solution g(£) = qo — 2d,t 2 + ■ ■ • . Henceforth we will 

consider only the large CI limit. The minimum of the potential at q m in is reached at time t m in = \/2d{qo — q m in)M — 
^(logQ)^ ~ * f° r l ar S e ^- Focussing then on the large time regime we can approximate the equation of motion by 

d 2 q {d-l) dq 

dt 2 + t dt ~ q 

Remarkably, for d = 3, we have an exact solution q(t) = f e _t with constant c fixed by matching to the small time 
solution so that q(t) = 5mm'min^"' f " tmi "' . Putting this solution into the action we find 

s (logClf 



log log CI 
for large CI. 

For d = 2, we mention for what it is worth that the equation of motion has an exact solution for large time in terms 
of a Bessel function: 

q(t) = c'K (t) ► dMe~\\ - 1 + • • • ) (4.14) 

t^oo y it at 

In fact, this instanton calculation seems to agree with the heuristic argument, given earlier in this section, only 
in d = 1. While we do not understand this discrepancy, we suspect that it is related to the existence of ultraviolet 



divergences for d > 1. Note that the heuristic estimate of Eq. (4.2) involves /(0), which is ultraviolet divergent for 
d > 1. On the other hand, the instanton result appears to be ultraviolet finite. However, this is somewhat illusory, 
since we expect loop corrections to bring in divergences, for d > 1. This issue requires further investigation. 

V. MAGNETIC FIELD 

In this paper, we propose a natural generalization of model I to include a magnetic field. In principle, this would 
define a new model for studying localization properties in the quantum Hall system. 
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There are two possibilities. In the first, we write 



f [Xi^Xjj — (J (.7*2 . Xj ) (i 



f !Cj A.dl 



(5.1) 



To have a well defined model, we have to specify the path joining Xi to Xj. The simplest choice is to take a straight 
line, then the phase factor in (5 ; l)becomes ei iB ^ XiX ^foi a co nsta nt magnetic field B. 



Unfortunately, f^ 1 {x,y), the functional inverse of f(x,y) in (5.1), does not have a particularly simple form, and so 
the corresponding field theory is not particularly attractive. 

On the other hand, we can simply define the model we would like to study by writing down the field theory 



S = 



» n 

j d 2 x J2(\D^ a \ 2 + m 2 \^ a \ 2 ) - pe-^:=^ 



(5.2) 



with the covariant derivative Dj = dj — iAj. Now the inverse f(x,y) of f~ x {x,y) — (-D 2 + m 2 )5(x — uLdoes not 
have a particularly simply form. The second form, Eq. (5.2) arises from the magnetic version of model II.E3 

Thus, we can define two different classes of models to study density of states and localization in the presence of a 
magnetic field. We can either have a simple f(x,y) or a simple f~ 1 (x,y). 

Suppose we want to calculate the density of states for p large. Doing a high density expansion of the Green's 
function, we have 



oo „ 

G ^ = E — ^ l / 



dxidx 2 :.dx n f(xi ~ xi)...f(xn-i - x n )e 



j B [(si X x 2 ) + {x2 X ) + . . . + (x~ n X si) ] 



After Fourier transforming 



fix) 



d 2 k 
(2tt)2 



7(fc) 



we can do the Gaussian integrals over x to obtain 



n=0 



d 2 k\ d 2 k n 



(2tt) 2 " (2tt) : 



f(h)...f (kn)e 



where 



^=E 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



where pj = kj — kj—t. We were not able to evaluate G(z) but in principle it might be possible for a Gaussian 
f(k) = e~ ak2 . 

In the B — > limit the stationary phase requirement forces all the fc's to be equal and so we recover the appropriate 
zero magnetic field result in Ref. ( |l|). In the opposite B — > oo limit we drop the exponential and obtain to leading 



order G(z) 



a result as expected. It simply says that as the Larmor radius goes to zero, each point in our 



z-pf(x=0) | 

random collection of points is its own universe. Again, in principle, the corrections in powers of X/B could be worked 
out. 
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